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Abstract. I briefly review our current understanding of dark matter and dark en- 
ergy. The first part of this paper focusses on issues pertaining to dark matter includ- 
ing observational evidence for its existence, current constraints and the 'abundance 
of substructure' and 'cuspy core' issues which arise in CDM. I also briefly describe 
MOND. The second part of this review focusses on dark energy. In this part I dis- 
cuss the significance of the cosmological constant problem which leads to a predicted 
value of the cosmological constant which is almost 10^^^ times larger than the ob- 
served value A/^-kG ~ 10~*^GeV'*. Setting A to this small value ensures that the 
acceleration of the universe is a fairly recent phenomenon giving rise to the 'cosmic 
coincidence' conundrum according to which we live during a special epoch when 
the density in matter and A are almost equal. Anthropic arguments are briefiy dis- 
cussed but more emphasis is placed upon dynamical dark energy models in which 
the equation of state is time dependent. These include Quintessence, Braneworld 
models, Chaplygin gas and Phantom energy. Model independent methods to deter- 
mine the cosmic equation of state and the Statefinder diagnostic are also discussed. 
The Statefinder has the attractive property a jali^ = 1 for LCDM, which is helpful 
for differentiating between LCDM and rival dark energy models. The review ends 
with a brief discussion of the fate of the universe in dark energy models. 
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1 Dark Matter 

Observations of the cosmic microwave background (CMB) and the deuterium 
abundance in the Universe suggest that J7baiyon/i^ — 0.02, or i7baryon — 0.04 
if the current Hubble expansion rate is h = iJo/lOOkm/sec/Mpc = 0.7. Al- 
though i7baryon is much larger than the observed mass in stars, i7stars — 0.005 
^, it is nevertheless very much smaller than the total energy density in the 
universe inferred from the observed anisotropy in the cosmic microwave back- 
ground |193| 

/2totai^^^g|^ = 1.02 ±0.02. (1) 

Both dark matter and dark energy are considered essential missing pieces in 
the cosmic jigsaw puzzle 

■^total '^baiyoiis — ^ (2) 

Although the nature of neither dark matter (DM) nor dark energy (DE) is 
currently known, it is felt that both DM and DE are non-baryonic in origin, 
and that DM is distinguished from DE by the fact that the former clusters on 
sub-Megaparsec scales (in order to explain galactic rotation curves) whereas 
the latter has a large negative pressure (and can make the universe accelerate) . 
In addition there is strong evidence to suggest that 

^ 1/3, r^DE ^ 2/3 . (3) 

In these lectures I will briefly review some properties of both dark matter and 
dark energy. 

Though the observational evidence favouring a flat Universe with i^totai — 
1 is fairly recent, the nature of the 'unseen' component of the universe (which 
dominates its mass density), is a long-standing issue in modern cosmology. 
Indeed, the need for dark matter was originally pointed out by Zwicky (1933) 
who realized that the velocities of individual galaxies located within the Coma 
cluster were quite large, and that this cluster would be gravitationally bound 
only if its total mass substantially exceeded the sum of the masses of its 
component galaxies. For clusters which have relaxed to dynamical equilibrium 
the mean kinetic and potential energies are related by the virial theorem [SHI 

K + ^ = 0, (4) 

where U ~ —GAP/R is the potential energy of a cluster of radius R, K ~ 
3M{vf)/2 is the kinetic energy and {v'^y^ is the dispersion in the line-of- 
sight velocity of cluster galaxies. (Clusters in the Abell catalogue typically 
have R ~ 1.5h^^ Mpc.) This relation allows us to infer the mean gravitational 

^ This suggests that most of the baryonic matter at 2: = is not contained in stars 
but might be contained in hot gas f30| . 
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potential energy if the kinetic energy is accurately known. The mass-to-light 
ratio in clusters can be as large as M/L ~ 300Mq/Lq. However since most 
of the mass in clusters is in the form of hot, x-ray emitting intracluster gas, 
the extent of dark matter in these objects is estimated to be M/Mium — 20, 
where Mium is the total mass in luminous matter including stars and gas. 



V (km/s) 




100 


* ■J-'^ expected 
- -...^ from 


. ■ 50 


r. ■ , luminous disk 




■ ■ 






M33 rotation curve 



Fig. 1. The observed rotation curve of the dwarf spiral galaxy M33 extends consid- 
erably beyond its optical image (shown superimposed); from Roy IIMI . 

In individual galaxies the presence of dark matter has been convincingly 
established through the use of Kepler's third law 

..(-) ^ ^ (5) 

to determine the 'rotation curve' v{r) at a given radial distance from the 
galactic center. Observations of galaxies taken at distances large enough for 
there to be no luminous galactic component indicate that, instead of declining 
at the expected rate v oc r~^/^ true if Af ~ constant, the velocity curves 
flattened out to u ~ constant implying M{r) (x r (see fignj. This observation 
suggests that the mass of galaxies continues to grow even when there is no 
luminous component to account for this increase. Velocity curves have been 
compiled for over 1000 spiral galaxies usually by measuring the 21 cm emission 
line from neutral hydrogen (HI) 148, ISTI. The results indicate that M/L = 
(10 — 20)Mq/Lq in spiral galaxies and in ellipticals, while this ratio can 
increase to M/L ~ (200 — 6OO)Mq/L0 in low surface brightness galaxies 
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(LSB's) and in dwarfs. For instance, a recent measurement of the Draco dwarf 
spheroidal galaxy located at a distance of only 79 kpc from the Milky Way 
shows the presence of a considerable amount of dark matter M/i| Draco = 
440 ± 240Mq /Lq ! It is interesting that the total mass of an individual 
galaxy is still somewhat of an unknown quantity since a turn around to the 
V oc r~^^'^ law at large radii has not been convincingly observed. 

An important difference between the distribution of dark matter in galax- 
ies and clusters needs to be emphasised: whereas dark matter appears to 
increase with distance in galaxies, in clusters exactly the reverse is true, the 
dark matter distribution actually decreases with distance. Indeed, for certain 
dwarfs (such as DD0154) the rotation curve has been measured to almost 15 
optical length scales indicating that the dark matter surrounding this object 
is extremely spread out (see also figure A foreground cluster, on the other 
hand, acts as a gravitational lens which focuses the light from background ob- 
jects such as galaxies and QSO's thereby allowing us to determine the depth of 
the cluster potential well. Observations of strong lensing by clusters indicate 
that dark matter is strongly concentrated in central regions with a projected 
mass of 10^'^ — 10^'* being contained within 0.2 - 0.3 Mpc of the central 
region. As we shall see later, this observation may prove to be problematic for 
alternatives to the dark matter hypothesis such as the Modified Newtonian 
Dynamics (MOND) approach of Milgrom [H^ . 

As discussed earlier, the fact that only 4% of the cosmic density is bary- 
onic suggests that the dark matter which we are observing could well be non- 
baryonic in origin. The need for non-baryonic forms of dark matter gets indi- 
rect support from the fact that baryonic models find it difficult to grow struc- 
ture from small initial conditions and hence to reconcile the existence of a well 
developed cosmic web of filaments, sheets and clusters at the present epoch 
with the exceedingly small amplitude of density perturbations {5p/p ~ 10~^ 
at z ~ 1, 100) inferred from COBE measurements and more recent CMB ex- 
periments. Indeed, it is well known that, if the effects of pressure are ignored, 
linearized density perturbations in a spatially flat matter dominated universe 
grow at the rate S oc i^/^ oc (1 -I- z)~^, where 1 + z — ao/a{t) is the cosmolog- 
ical redshift. (Contrast this relatively slow growth rate with the exponential 
'Jeans instability' of a static matter distribution S oc exp y/^nGpt.) In a bary- 
onic universe, the large radiation pressure (caused by thompson scattering of 
CMB photons off electrons) ensures that density perturbations in the baryonic 
component can begin growing only after hydrogen recombines at z ~ 1, 100 
at which point of time baryons and radiation decouple. Requiring 5 > 1 today 
implies 5 > 10~^ at recombination, which contradicts CMB observations by 
over an order of magnitude ! In non-baryonic models on the other hand, the 
absence of any significant coupling between dark matter and radiation allows 
structure to grow much earlier, significantly before hydrogen in the universe 
has recombined. After recombination baryons simply fall into the potential 
wells created for them by the dominant non-baryonic component. As a result 
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a universe with a substantial non-baryonic component can give rise to the 
structure which we see today from smaller initial fluctuations. 

The growth of structure via gravitational instability depends both upon 
the nature of primordial perturbations (adiabatic/isocurvature) and upon 
whether the dark matter species is hot or cold. The issue of density per- 
turbations has been discussed in considerable detail by Ruth Durrer at this 
school and I shall not touch upon this important topic any further. Let me 
instead say a few words about hot and cold dark matter. Non-baryonic Hot 
Dark Matter (HDM) particles are assumed to have decoupled from the rest of 
matter/radiation when they were relativistic and so have a very large veloc- 
ity dispersion (hence called 'hot'). Cold Dark Matter (CDM) particles, on the 
other hand, have a very small velocity dispersion and decouple from the rest of 
matter/radiation when they are non- relativistic. The free-streaming (collision- 
less phase mixing) of non-baryonic particles as they travel from high density 
to low density regions (and vice versa) introduces an important length scale 
called the 'free-streaming distance' Afs - which is the mean distance travelled 
by a relativistic particle species until its momentum becomes non-relativistic. 
In both HDM and CDM the processed final spectrum of density perturbations 
differs from its initial form. In the case of HDM this difference arises because 
fluctuations on scales smaller than Afs are wiped out due to free streaming 
with the result that the processed final spectrum has a well defined cutoff on 
scales smaller than A ^ Afs. Perhaps the best example of HDM is provided 
by a hght neutrino of mass about 30 eV. In this case Afs ~ 41(30eV/mi,) 
Mpc with the result that large proto-pancakes having masses comparable to 
those of rich clusters of galaxies M ^ IO^^Mq are the first objects to form 
in HDM. Smaller objects (galaxies) are formed by the fragmentation of the 
proto-pancake. This top-down scenario for structure formation was originally 
suggested by Zeldovich and coworkers in connection with adiabatic baryonic 
models and subsequently applied to HDM. It has since fallen out of favour 
mainly due to the strong observational constraints on the mass of the neu- 
trino J2iy < 0.7 eV and on the rehc neutrino density 10"'^ < H^h'^ < 10^^ 
| 61l 11931 1631 1123| . It also faces considerable difficulty in forming structure suf- 
ficiently early to explain the existence of galaxies and QSO's at high redshifts. 

In contrast to HDM, constituents of CDM have a much smaller free- 
streaming distance. Because of this small scales are the first to go non-linear 
and gravitational clustering proceeds in a bottom up fashion in this scenario. 
A key quantity defining gravitational clustering is the power spectrum of den- 
sity perturbations P{k) = \Sk\'^, which is related to the mean square density 
fiuctuation via 



Inflationary models predict Pi(fc) cx fe", ri ~ 1, at an early epoch. As the uni- 
verse expands the power spectrum gets modified. The 'processed' final spec- 
trum depends upon the nature of dark matter, the epoch of matter-radiation 




(6) 
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equality and other cosmological quantities. The final and initial spectra are 
related through a transfer function 

Pf{k)^P,{k)xT^k) . (7) 

CDM-type spectra have the following approximate form of the transfer func- 
tion [TCSIITMIITC^ 

Equations Q & (|HJ| illustrate the 'turn around' of the power spectrum from 
its primordial scale invariant form P{k) cx fc on the largest scales to P{k) oc 
k~^ log^ k on small scales. (The precise location of the turn-around and the 
amplitude of P{k) depend upon specific details of the cosmological model, see 
for instance 

The 'standard' cold dark matter (SCDM) paradigm, which assumed that 
■i^CDM — 1, was introduced during the early 1980's at roughly the same time 
when HDM was perceived to be in trouble (see |l(Jlll96[|5nill66| for references 
to earlier work on this subject). Although SCDM was very successful in ex- 
plaining a host of observational details, it was clear already a decade ago, that 
the processed power spectrum of SCDM lacked sufficient power on large scales 
and so fell short of explaining the angular two point correlation function for 
galaxies on scales ^ 50 Mpc [201. The relevant cosmological quantity in this 
case is the shape of the power spectrum of density perturbations, which for 
CDM-like models, can be characterised by the 'shape parameter' F — Hj^h. 
SCDM models with — 1 and the HST-determined value h ~ 0.7 predict 
r ~ 0.5 which is much larger than the observed value P = 0.207 ± 0.030 
inferred from observations of galaxy clustering in the sloan digital sky survey 
(SDSS) 154 . A modification of SCDM called LCDM assumes that, in addi- 
tion to CDM the universe consists of a smoothly distributed component called 
a cosmological constant or a Lambda-term. LCDM models with h ~ 0.7 and 
~ 0.3 predict a smaller value for the shape parameter, P ~ 0.2, and the 
resulting amplitude and shape of the power spectrum is in excellent agreement 
with several different sets of observations as demonstrated in figure [3 

From ijBJ, ITJ & (jSJ we find that on small scales, the contribution to the 
rms density fluctuation from a given logarithmic interval in k is 

(^^)^^fc3p^(fc)(xlog^A , (9) 

which illustrates the fact that, although the smallest scales are the first to go 
non-linear, there is significant power to drive gravitational instability rapidly 
to larger scales in this model. Indeed, detailed N-body simulations of large 
scale structure show that filaments defining the cosmic web first form on the 
smallest scales. The scale-length characterizing the cosmic web grows as the 
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Fig. 2. The power spectrum inferred from observations of large scale structure, 
the LymauQ forest, gravitational lensing and the CMB. The solid line shows the 
power spectrum prediction for a flat scale-invariant LCDM model with J7m = 0.28, 
QilQm = 0.16, h = 0.72; from Tegmark et al. {200^ . 

universe expands, until at the present epoch the cosmic web consists of a 
fully developed supcrclustcr-void network with a scale-length of several tens 
of Megaparsec ^81, 183, 119, 206 . 

Promising candidates for cold dark matter include a 100— 1000 GeV parti- 
cle called a neutralino. The neutralino is a weakly interacting massive particle 
(WIMP). As its name suggests it is neutral and is a fermionic partner to 
the gauge and Higgs bosons (usually called the 'bino, wino and higgsino'). It 
is believed that the lightest supersymmetric particle will be stable due to R- 
parity which makes the neutralino an excellent candidate for cold dark matter 
(see [ 1631 189| for reviews of particle dark matter). A radically different particle 
candidate for cold dark matter is an ultra-light pseudo-Goldstone boson called 
an axion with a mass of only ~ 10~^^^ eV. Although ultralight, the axion 
is 'cold' because it was created as a zero-momentum condensate. Its existence 
is a by-product of an attempt to resolve QCD of what is commonly called the 
'strong CP problem' which arises because non-perturbative effects in QCD 
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give rise to an electric dipole moraent for the neutron - in marked contrast 
with observations |1(J1| . Other candidates for non-baryonic cold dark matter 
include string theory motivated modulii fields |32j : non-thermally produced 
super-heavy particles having a mass ~ 10^'' GeV and dubbed Wimpzillas [TnU| : 
as well as axino's and gravitino's - superpartners of the axion and graviton 
respectively \l(y'6\ . 

Since WIMP's cluster gravitationally, one should expect to find a flux of 
these particles in our own solar system and attempts are being made to deter- 
mine dark matter particles by measuring the scattering of WIMP's on target 
nucleii through nuclear recoils. Now the earth orbits the sun with a velocity 
~ 30 km/sec, even as the sun orbits the galaxy with vi\i^^ ~ 220 km/sec. 
Furthermore the plane of the Earth's orbit is inclined at an angle of 60° to 
the glactic plane, because of which the dark matter flux on Earth is expected 
to be larger in June (when the Earth's velocity and the Sun's velocity add 
together) than in December (when these two velocities subtract). The result- 
ing rate variation is about 7% between the flux measured during summer and 
winter. Precisely such a signal was reported by the DAMA experiment whose 
data (collected since 1996) appears to show a yearly modulation with greater 
events reported in June than in December |17| . However results obtained by 
the DAMA group remain controversial since they have not been substanti- 
ated by other groups which report negative results for similar searches (see 
fl29.„.95j for recent reviews on this subject). 




Fig. 3. The Earth's motion around the Sun; from Klialil and Munoz (2001). 



Despite the excellent agreement of LCDM with large scale observations, 
some concerns have recently been expressed about the ability of this model to 
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account for a number of smaller scale observations which can be summarized 
as follows: 

• The substructure problem is used to describe the fact that the cold dark 
matter model (with or without a cosmological constant) predicts an ex- 
cessive number of dark matter subhaloes (or substructure) within a larger 
halo. If one (perhaps naively) associates each halo with a gravitation- 
ally bound baryonic object then the predicted number of dwarf-galaxy 
satellites within the local group exceeds the observed number by over 
an order of magnitude. Indeed, detailed N-body simulations as well as 
theoretical estimates predict around 1000 dark matter satellites in our 
local group which is much larger than the 40 or so observed at present 

[na iiii El 133 1112 El inni iins Ei ■ 

• The cuspy core problem 

CDM predicts a universal density profile for dark matter halos in the 
wide range 1O^M0 — IO^^Mq which applies both to galaxy clusters as well 
as individual galaxies including dwarfs and LSB's. ^ The density profile 
originally suggested by Navarro, Frenk and White |133j is 

(10) 

which gives p oc for r <C and p oc r^^ for r 3> r^, where is the 
scale radius and po is the characteristic halo density. (Other groups using 
higher resolution computations found somewhat steeper density profiles at 
small radii, such as p oc r~^-^ |127[I57] .) 

The cuspy core problem refers to the apparent contradiction between N- 
body experiments - which show that the density profile in CDM halos has 
a 1/r (or steeper) density cusp at the center, and observations - which 
appear to favour significantly shallower density cores in galaxy clusters as 
well as in individual dwarf and LSB galaxies (see 1^ 1571 [T^ [T^ 

UnmHHnnHinil for detailed discussions of this issue). 

Although disconcerting, given the very considerable success of LCDM in 
explaining gravitational clustering on large scales, it may at this point be pre- 
mature to condemn this model on the basis of small scale observations alone. 
It could be that the difRculties alluded to above are a result of an oversim- 
plification of the complex physical processes involved and that a more careful 
analysis of the baryonic physics on small scales including the hydrodynamical 
effects of star formation and supernova feedback needs to be undertaken. For 
instance both dwarfs and LSB's have very shallow potential wells, a strong 
burst of star formation and supernova activity may therefore empty dark 

^ Low Surface Brightness Galaxies (LSB's) are dominated by their dark matter 
content and therefore provide particularly good astrophysical objects with which 
to test dark matter models. 



p{r) = po [rjr) 



l+(^ 
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matter halos of their baryonic content resulting in a large number of 'failed 
galaxies' and providing a possible resolution to the 'satellite catastrophe'. 
(The failed galaxies will act as gravitational lenses and should therefore be 
detectable through careful observations.) Other explanations include the ef- 
fects of tidal stripping recently discussed in |103j . Likewise issues involving 
beam smearing, the influence of bars and the interaction of baryons and dark 
matter in the central regions of galaxies and clusters could be intricately linked 
with the central cusp issue and must be better understood if one wishes to 
seriously test the CDM hypothesis on small scales. 

In concluding this discussion on dark matter I would like to briefly mention 
Modified Newtonian Dynamics (MOND) which, in some circles, is regarded as 
an alternative to the dark matter hypothesis. As the name suggests, MOND is 
a modification of Newtonian physics which proposes to explain the flat rota- 
tion curves of galaxies without invoking any assumptions about dark matter. 
Briefly, MOND assumes that Newtons law of inertia {F = ma) is modified at 
sufficiently low accelerations (a < uq) to 

F = ma^{a/ao) , (11) 

where ^(x) = x when x and ijl{x) — 1 when a; ^ 1 |122II18()] . It is easy to 
see that this results in the modification of the conventional formula for gravi- 
tational acceleration F — mgN, resulting in the following relation between the 
true acceleration and the Newtonian value: a = ^JgWo-o- For a body orbiting 
a point mass M, = GMjr^. Since the centripetal acceleration a = jr 
now equals the true acceleration a, one gets 

w'^ = GA/ao , (12) 

i. e. for sufficiently low values of the acceleration the rotation curve of an iso- 
lated body of mass M does not depend upon the radial distance r at which the 
velocity is measured, in other words not only does this theory predict flat rota- 
tion curves it also suggests that the individual halo associated with a galaxy 
is infinite in extent ! (This latter prediction may be a problem for MOND 
since recent galaxy-galaxy lensing results tH21 suggest that galaxy halo's may 
have a maximum extent of about 0.5 Mpc.) The value of ao needed to explain 
observations is ap ^ lO^^cm/s^ which is of the same order as ci?o ! This has 
led supporters of this hypothesis to conjecture that MOND may reflect "the 
effect of cosmology on local particle dynamics" ,180, . Although MOND gives 
results which are in good agreement with observations of individual galaxies, 
it is not clear whether it is as successful for explaining clusters for which strong 
gravitational lensing indicates a larger mass concentration at cluster centers 
than accounted for by MOND (THDIESI- Another difficulty with MOND is that 
it is problematic to embed this theory within a more comprehensive relativis- 
tic theory of gravity and hence, at present, it is not clear what predictions 
a MOND-type theory may make for gravitational lensing and other curved 
space-time effects. For some recent developments in this direction see |23| . 



Dark Matter and Dark Energy 



11 



To summarise, current observations make a strong case for clustered, non- 
baryonic dark matter to account for as much as a third of the total matter 
density in the Universe f2,„ ~ 1/3. The remaining two-thirds is thought to 
reside in a relative smooth component having large negative pressure and 
called Dark Energy. 



2 Dark Energy 

2.1 The cosmological constant and Vacuum energy 

Type la supernovae, when treated as standardized candles, suggest that the 
expansion of the universe is speeding up rather than slowing down. The case 
for an accelerating universe also receives independent support from CMB and 
large scale structure studies. All three data sets can be simultaneously satified 
if one postulates that the dominant component of the universe is relatively 
smooth, has a large negative pressure and I2de — 2/3. 

The simplest example of dark energy is a cosmological constant, intro- 
duced by Einstein in 1917. The Einstein equations, in the presence of the 
cosmological constant, aquire the form 

1 SttG 

^ik — -^gikR — ^4 Tik + Agik ■ (13) 

Although Einstein originally introduced the cosmological constant {A) into 
the left hand side of his field equations, it has now become conventional to 
move the yl-term to the RHS, treating it as an effective form of matter. In 
a homogeneous and isotropic Friedmann-Robertson- Walker (FRW) universe 
consisting of pressureless dust (dark matter) and A, the Raychaudhury equa- 
tion, which follows from 1)13(1 . takes the form 

47rG A 
a = ^apm + ■ (14) 

Equation p4|l can be rewritten in the form of a force law: 

+ (15) 

which demonstrates that the cosmological constant gives rise to a repulsive 
force whose value increases with distance. The repulsive nature of A could be 
responsible for the acceleration of the universe as demonstrated in ((14(1 . 

Although introduced into physics in 1917, the physical basis for a cosmo- 
logical constant remained a bit of a mystery until the 1960's, when it was 
realised that zero-point vacuum fluctuations must respect Lorenz invariance 
and therefore have the form (T^fc) = Agi]^ |214( . As it turns out, the vacuum 
expectation value of the energy momentum is divergent both for bosonic and 
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fermionic fields, and this gives rise to what is known as 'the cosmological 
constant problem'. Indeed the effective cosmological constant generated by 
vacuum fluctuations is 

^ = ('roo)vac « / Vk^ + m^k^dk , (16) 
SttG Jq 

since the integral diverges as fc^ one gets an infinite value for the vacuum 
energy. Even if one chooses to 'regularize' (Tik) by imposing an ultraviolet 
cutoff at the Planck scale, one is still left with an enormously large value for 
the vacuum energy (Too)vac — /G^H lO^^GeV* which is 123 orders of 
magnitude larger than the currently observed pA — 10^*^GeV"'. A smaller 
ultraviolet cut-off does not fare much better since a cutoff at the QCD scale 
results in Aq^^jj ^ 10~'^GeV'*, which is still forty orders of magnitude larger 
than observed. 

In the 1970's the discovery of supersymmetry led to the hope that, since 
bosons and fermions (of identical mass) contribute equally but with opposite 
sign to the vacuum expectation value of physical quantities, the cosmological 
constant problem may be resolved by a judicious balance between bosons 
and fermions in nature. However supersymmetry (if it exists) is broken at 
the low temperatures prevailing in the universe today and on this account 
one should expect the cosmological constant to vanish in the early universe, 
but to reappear during late times when the temperature has dropped below 
TsuSY- This is clearly an undesirable scenario and almost the very opposite 
of what one is looking for, since, a large value of A at an early time is useful 
from the viewpoint of inflation, whereas a very small current value of A is in 
agreement with observations [1721 11711 . 

In the absence of a resolution to the cosmological constant problem the 
following possibility connecting the vacuum energy with the SUSY and Planck 
scales may be worth exploring |172l I171j . The mass scale associated with 
the scale of supersymmetry breaking in some models, AfsuSY ~ 1 TeV, lies 
midway between the Planck scale and 10~^ eV. One could conjecture that 
the small observed value of the cosmological constant pA — (lO'^eV)"^ is 
associated with the vacuum in a theory which had a fundamental mass scale 
Mx Ml^sy/Mph such that pvac ~ - (IQ-^eV)^. 

The cosmological constant is also relevant from the perspective of mod- 
els with spontaneous symmetry breaking |2U9| . Indeed, if one examines the 
Lagrangian 

£ = i 5'-'9,</)a,0-y(0), 

F(0) ^Vo- + ^Xc^\ (17) 

one notices that the symmetric state at = is unstable and the system 
settles in the ground state </> = +(t or = — cr, where a = y'Ji^/X, thereby 
breaking the reflection symmetry <-> — present in the Lagrangian. For 
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Vq = this potential gives rise to a large negative cosmological constant 
Af,ff = V{(f> — a) ^ — /x*/4A in the broken symmetry state. This embarassing 
situation can be avoided only if one chooses a value for Vq which almost 
exactly cancels Ae/f, namely Vq — +/^''/4A so that Aeff/SnG = Vq — /x^/4A ~ 
10-47GeV 4. 

The cosmological consequences of this rather ad-hoc 'regularization' ex- 
ercise are instructive. Unless the value of Aeff lies in a very small window, 
the universe will be a very different place from the one we are used to. For 
instance a negative value of the yl-term yle/j/STrG < — 10~^'^GeV will cause 
the universe to recoUapse (the effect of A is attractive now instead of being 
repulsive) less than a billion years after the big bang - a period which is much 
too short for galaxies to form and for life (as we know it) to emerge. On the 
other hand a large positive yle///87rG > 10~'''^GeV'' makes the universe accel- 
erate much before the present epoch, thereby inhibiting structure formation 
and precluding the emergence of life. 

The very small window in A which allows life to emerge has led some 
cosmologists to propose anthropic arguments for the existence of a small cos- 
mological constant [23 IUSI 171)1 121()j . One such possibility is the following "if 
our big bang is just one of many big bangs, with a wide range of vacuum en- 
ergies, then it is natural that some of these big bangs should have a vacuum 
energy in the narrow range where galaxies can form, and of course it is just 
these big bangs in which there could be astronomers and physicists wondering 
about the vacuum energy" j21()| . 

I will not discuss the anthropic argument any further in these lectures but 
will point the interested reader to 118. I76[[21U| for further discussion of this 
issue. 

It is important to note that there is no known fundamental symmetry in 
nature which will set the value of A to zero. In its absence, the small observed 
value of the dark energy remains somewhat of a dilemma which remains to 
be fully understood and resolved. ^ 



2.2 Dynamical models of Dark Energy 

The cosmological constant is but one example of a form of matter (dark en- 
ergy) which could drive an accelerated phase in the history of our universe. 
Indeed, (fTHl is easily generalised to 

- ^Y.^P^ + 3p,) = —Y, p,{l + Zw,) , (18) 

i i 

where the summation is over all forms of matter present in the universe with 
equation of state Wi = Pi/ pi- Eqn. (|18(l together with its companion equation 

The important role played by symmetries is illustrated by the U(l) gauge symme- 
try of electrodynamics whose presence implies a zero rest mass for the photon. No 
analogous symmetry exists for the neutrino and recent experiments do indicate 
that neutrino's could have a small mass. 
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A < 



Fig. 4. Spontaneous symetry breaking in many field theory models takes the form 
of the Mexican top hat potential shown above. The dashed line shows the potential 
before the cosmological constant has been 'renormalized' and the solid line after. 
(From Sahni and Starobinsky 2000.) 



completely describes the dynamics of a FRW universe {k/a? is the Gaussian 
curvature of space). 

Clearly a universe consisting of only a single component will accelerate if 
w < —1/3. Fluids satisfying p + 3p > oi w > —1/3 are said to satisfy the 
'strong energy condition' (SEC). We therefore find that, in order to acceler- 
ate, 'dark energy' must violate the SEC. Another condition which is usually 
assumed to be sacrosanct, but has recently been called into question is the 
'weak energy condition' (WEC) p + p>Ooiw> —1. Failure to satisfy the 
WEC can result in faster-than-exponential expansion for the universe and in 
a cosmic 'Big Rip', which we shall come to in a moment. 

It is often more convenient to rewrite (|18|l in terms of the 'deceleration 
parameter' 




(19) 



q = 



am 



a 



3H^ 



(1 + 3wx^?x) 
2 



(20) 



where Hi — SnGpi /SH"^ and we have assumed a flat universe with [2^ + i^x — 
1 {Six = ^^de)- The condition for accelerated expansion {q < 0) is equivalent 
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to 



which leads to 



w < for n„,^Q , (22) 

w < -i for a,n = 1/3 . (23) 

Eqn. (|19|) can be used to develop an expression for the Hubble parameter 
H = a/a in terms of the cosmological redshift z = ao/a(i) — 1: 

H{z) = Ho + zf + nx{l + , (24) 

where Hq = H{z = 0) is the present value of the Hubble parameter, — 
SirGpom/SHQ, fix — 87rGpoDE/3i?0' describe the dimensionless density of 
matter and dark energy respectively, {w = wde)j and we have made the 
assumption of a flat universe so that ^?,„ + f2x = 1- 

In LCDM cosmology w — —1, SI a — A/3Hq, and the expansion factor has 
the elegant form ^172, 

a(0(x fsinh I ^|ct J , (25) 

which smoothly interpolates between a matter dominated universe in the past 
(a oc i^/^) and accelerated expansion in the future (a cx exp \J^t). 

We are now in a position to appreciate the evidence for an accelerating 
universe which originates in observations of the light flux from high redshift 
type la supernovae. Type la supernovae are extremely bright objects, {Mb — 
— 19.5) which makes them ideally suited for studying the properties of the 
universe at large distances. 

The light flux received from a distant supernova is related to its absolute 
luminosity C and its 'luminosity distance' cIl through the relation 

If one views this problem from within the Newtonian perspective then, since 
the geometry of space is Euclidean, c^l = ^/x^ + y^ + z^. In general relativity, 
on the other hand, the geometry of space can be non-Euclidean, and the 
luminosity distance to an object located at redshift z will, in general, depend 
both upon the geometry of space as well as the expansion history of the 
universe. Indeed, it can be shown that in a spatially flat and expanding FRW 
universe, the luminosity distance has the form 
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ddz)-il + ^)l J^y (27) 

The luminosity distance is shown in Fig. 1 for a number of cosmological 
models with varying amounts of f2m & ^A- The limiting case f^m = 1, = 
corresponds to standard cold dark matter (SCDM) in which the universe 
decelerates as a weak power law a{t) oc t"^/^. The other extreme example 

= 1, Qm = describes the de Sitter universe (also known as steady state 

cosmology) which accelerates at the steady rate a{t) cx exp ij^t. From figure 
1 we see that a supernova at redshift z = 3 will appear 9 times brighter in 
SCDM than it will in de Sitter space ! 



o 




Z 4 6 8 10 

redshift z 



Fig. 5. The luminosity distance dL (in units of Hq^) is shown as a function of 
cosmological redshift z for spatially flat cosmological models with J7m + Oa = 1- 
Heavier lines correspond to larger values of ilm- The dashed line shows the lumi- 
nosity distance in the spatially flat de Sitter universe {SI a = 1). From Sahni and 
Starobinsky (T7^ . 



Systematic studies of type la supernovae have revealed that: 

• Type la Sn are excellent stadardized candles. The dispersion in peak su- 
pernova luminosity is small: Am ~ 0.3, and the corresponding change in 
intensity is about 25%. In addition the light curve of a type la supernova 
is correlated with its peak luminosity 149 to a precision of ~ 7%, so 
that brighter supernovae take longer to fade. (Type la Sn take roughly 
20 days to rise from relative obscurity to maximum light.) This allows 
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us to 'standardize' supernova light curves thereby reducing the scatter in 
their luminosities to 12% which turns type la supernovae into very good 
standard candles. 

• Type la supernovae at higher redshifts are consistently dimmer than their 
counterparts at lower redshifts relative to what might be expected in 
SCDM cosmology. If type la supernovae are treated as standard candles 
then, assuming systematic effects such as cosmological evolution and dim- 
ming by intergalactic dust are either not vitally important or have been 
corrected for, the systematic dimming of high— z Sn can be interpreted as 
evidence for an accelerated expansion of the universe caused by a form of 
'dark energy' having large negative pressure. 

The evidence for an accelerating universe from high redshift type la super- 
novae has now received independent support from an analysis of CMB fluc- 
tuations together with the HST key project determination of the Hubble pa- 
rameter. Interestingly, the degeneracy in parameter space {f2rm ^a} arising 
from Sn observations is almost orthogonal to the degeneracy which arises from 
CMB measurements. This principle of 'cosmic complementarity' serves to sig- 
nificantly reduce the errors on J7„j & 12^1 when the two sets of observations 
are combined, as shown in figure El 




Fig. 6. Constraints on the density of dark matter Om and dark energy in the form of 
a cosmological constant Oa, determined using WMAP (upper left), WMAP -I- other 
CMB experiments (WMAPext; upper right), WMAPext + HST key project data 
(lower left) and WMAPext -I- HST -I- supernova data (lower right); from Spergel et 
al (2003). 
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If dark energy is described by an unevolving equation of state w ~ px / Px i 
then the transition between deceleration and acceleration (a = 0) occurs at 
the redshift 

{l + z^y^^ = -{l + 3w)^ w <0. (28) 

Another important redshift describes the epoch when the densities in dark 
matter and dark energy are equal 

(l + -oq)^'"-(^) . (29) 

Substituting Qa = 0.7, Qm = 0.3 we find ~ 0.73, z^q ~ 0.37 for LCDM. 
The fact that the acceleration of the universe is a fairly recent phenomenon 
illustrates the 'cosmic coincidence' puzzle according to which we appear to live 
during a special epoch when the densities in dark energy and in dark matter 
are almost equal. A recent origin for the acceleration epoch is supported by 
supernova observations which suggest a decelerating universe aX z >^ 0.5 |16U| . 
It is important to note that dark enegy models with an unevolving equation of 
state need to have their initial conditions properly 'tuned' in order to dominate 
the universe at precisely the present epoch. This problem is most acute for the 
cosmological constant. Since the cosmological constant does not evolve while 
both matter and radiation evolve rapidly {pm oc a~^, Pr oc a"'*), it follows that 
the small current value pA = A/^nG ~ 10~^^ GeV* implies Pa/ Pr — 10"^^'^ 
at the Planck time (when the temperature of the universe was T ~ 10^^ GeV), 
or paI Pr — 10~^^ at the time of the electroweak phase transition (T ~ 100 
GeV). Thus an extreme fine-tuning of initial conditions is required in order 
to ensure that pA / pm ^ 1 today ! 

The fine tuning problem which plagues A also affects DE models in which 
w = constant ^ — 1. A combined analaysis of CMB, galaxy clustering and 
supernovae data indicates that a constant equation of state for dark energy 
must satisfy w < —0.82 at the 95% confidence level |193l I2U1) . and it is 
easy to show that for these models the fine tuning (and cosmic coincidence) 
problems are almost as acute as they are for the cosmological constant. This 
constraint on w also virtually rules out two interesting DE candidates based 
on topological defect models: a tangled network of cosmic strings w ~ —1/3 
and domain walls w ~ —2/3. 

2.3 Quintessence 

It is interesting that the fine tuning problem facing dark energy models with 
a constant equation of state can be alleviated if we assume that the equation 
of state is time dependent. An important class of models having this property 
are scalar fields (quintessence)^ which couple minimally to gravity so that 

* Quintessence is named after the all pervasive fifth element of ancient philosophical 
thought. Note that the quintessence Lagrangian is the same as that used for 
Inflationary model building. 
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their Lagrangian density and energy momentum tensor is 

£ = _ ^(^) ^ (30) 

P^TS = 1<P' + V{cb), p ^ -T^ = - V{<l^) , (31) 

where we have assumed, for simpUcity, that the field is homogeneous. Po- 
tentials which are sufficiently steep to satisfy 7^ = V"V / {VY > 1 have the 
interesting property that scalar fields rolling down such a potential approach 
a common evolutionary path from a wide range of initial conditions |217j . In 
these so-called 'tracker' models the scalar field density (and its equation of 
state) remains close to that of the dominant background matter during most 
of cosmological evolution. 




z+1 

Fig. 7. The quintessence Q-field while rolling an inverse power law potential tracks 
first radiation then matter, before coming to dominate the energy density of the 
universe at present. If the initial value of the Q-field density is small then pq remains 
constant until pq ~ Prad, and then follows the tracker trajectory. From Zlatev, Wang 
and Steinhardt PT7| . 



An excellent example of a tracker potential is provided by V{(p) = Vo/(/)" 
| 157| . During tracking the ratio of the energy density of the scalar field 
(quintessence) to that of radiation/matter gradually increases p^/pB c<. 
^4/{2+a) i^g equation of state remains marginally smaller than the back- 

ground value = {awB — 2)/{a + 2). For large values of (j) this potential 
becomes flat ensuring that the scalar field rolls sufficiently slowly (0^ ^ y{4>)) 
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to allow the universe to accelerate. Note that for quintessence fields the con- 
dition's 1221) & translate into 

w<p<-^ ^ 0^ < 1^(0) , 

W4.<-^ =^ 0' < lv{4>) . (32) 

(Current observations imply a < 2.) 

An extreme example of quintessence is provided by the exponential poten- 
tial V{(f>) = Vo exp (- V8^A0/Mp) [1231213, where Mp = 1/VG is the Planck 
mass. In this case 

— = -5 = constant < 0.2, (33) 

Pb + P0 

is the background energy density while wb is the associated background 
equation of state. The lower limit /C^/ptotai < 0.2 arises because of nucleosyn- 
thesis constraints which prevent the energy density in quintessence from being 
large initially (at t ^ few sec). Equation (|33|l suggests that the exponential 
potential will remain subdominant if it was so initially. An interesting poten- 
tial which interpolates between an exponential and a power law can however 
give rise to late time acceleration from tracker-like initial conditions |168| 

V{(j)) = Vo [cosh A</> - 1]P, (34) 

has the property that ~ wb at early times whereas {w^) — (p — l)/(p + 
1) at late times. Consequently ()34|l describes quintessence for p < 1/2 and 
pressureless 'cold' dark matter (CDM) for p = 1. Thus the cosine hyperbolic 
potential (|34|l is able to describe both dark matter and dark energy within a 
tracker framework (also see |204l I12| '). 

Remarkably, quintessence can even accomodate a constant equation of 
state {w = constant) by means of the potential (1721 11731 12l)3j 

V{(j)) cx sinh^^ {C4> + D) , (35) 
with suitably chosen values of C, D. 

Quintessence models can be divided into two categories: models which 
roll to large values of (f>/mp >, 1 and models for which (f>/mp <C 1 at the 
present epoch. An important concern for the former is the effect of quantum 
corrections which, if large, could alter the shape of the quintessence potential 
[ins ESI EHI UBS- An important related issue is that the coupling between 
standard model fields and quintessence must be small in order to have evaded 
detection. Moreover even small couplings between quintessence and standard 
model fields can give rise to interesting changes in cosmology as shown in 
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Quintessence Potential Reference 



Voexp(-A0) Ratra & Peebles (1988), Wetterich (1988), 

Ferreira & Joyce (1998) 

m^0^, A(^* Frieman et al (1995) 

Vo/^", a > Ratra & Peebles (1988) 

Vb exp (At/)^)/^" Brax & Martin (1999,2000) 

V^)(cosli A0 - 1)P Sahni & Wang (2000) 

Vosinh"" (A0) Sahni & Starobinsky (2000), Urefia-Lopez & Matos (2000) 

Vo{e°""'' + e^'"'') Barreiro, Copeland & Nunes ( 2000) 

V'o(exp Mp/<^ - 1) Zlatev, Wang & Steinhardt (1999) 

yo[(0 - B)" + AJe"^"* Albrecht & Skordis (2000) 

Table 1. 



I would like to end this section by mentioning that, due to the shortage of 
time I have not been able to cover all of the DE models suggested in the liter- 
ature (a number that is growing rapidly !) For this reason these lectures will 
not discuss DE due to vacuum Dolarization |167[ [141] . k-essence ITS', Cardas- 
sian expansion [73j, Quasi-Steady State Cosmology |132| . scalar-tensor models 
IP", "^S", U^, TM IMl 1^ 1113 • For other interesting approaches see 
[TH1IE3, _8i 12 , 137, SQiUnmillinilllHZIEISl- a partial list of some popular 
quintessence models is given in Table 1, and the reader is also referred to the 
dark energy reviews in [TTl HTl [1151 ITTTl [HOI . 



2.4 Dark energy in braneworld models 

Inspired by the Randall- Sundrum |156j scenario, braneworld cosmology sug- 
gests that we could be living on a three dimensional 'brane' which is embedded 
in a higher (usually four) dimensional bulk. According to such a scheme, all 
matter fields are confined to the brane whereas the graviton if free to pro- 
pogate in the brane as well as in the bulk (see the lectures by Roy Maartens 
in this volume and |110| for a comprehensive discussion of Braneworld cos- 
mology.) Within the RS setting the equation of motion of a scalar field pro- 
pogating on the brane is 
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4> + 3H(P + V'{(j)) ^0. (36) 



where [TM] 



3m^ 2(7 3 

P = ^-^^ + V{4>) ■ (37) 

f is an integration constant which transmits bulk graviton influence onto the 
brane. The brane tension a provides a relationship between the four dimen- 
sional Planck mass (m) and the five-dimensional Planck mass (M) 



(38) 

a also relates the four-dimensional cosmological constant A4 on the brane to 
the five-dimensional (bulk) cosmological constant Ah through 

47r / , 47r ,\ , , 

^4 = T7?Mb + x-F^^M- (39) 




A'py " 3M3 

Note that H37|l contains an additional term p^/cr whose presence can be 
attributed to junction conditions imposed at the bulk-brane boundary. Be- 
cause of this term the damping experienced by the scalar field as it rolls down 
its potential dramatically increases so that infiation can be sourced by poten- 
tials which are normally too steep to produce slow-roll. Indeed the slow-roll 
parameters in braneworld models (for V/cr ^ 1) are [ITT] 

e ~ 4eFRw( W\ V - 2j^Fnw{V/a)-\ (40) 

illustrating that slow- roll (e,7y <C 1) is easier to achieve when V/cr 3> 1. 
Infiation can therefore arise for the very steep potentials associated with 
quintessence such as F (x e~^'^, V oc etc. This gives rise to the in- 

triguing possibility that both infiation and quintessence may be sourced by 
one and the same scalar field. Termed 'quintessential infiation', these models 
have been examined in [HI ESI ESI HHi 1111 EHl 11001 E3 CHI CZHl IHTj. An 
example of quintessential inflation is shown in figure El 

A radically different way of making the Universe accelerate was sug- 
gested in |55[ I170j . The braneworld model developed by Deffayet, Dvali and 
Gabadadze (DDG) was radically different from the RS model in that both 
the bulk cosmological constant and the brane tension were set to zero, while 
a curvature term was introduced in the brane action so that the theory was 
described by 

S^M^ [ n + m" I R + I L„,atter ■ (41) 
J bulk J branc J brane 
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Fig. 8. The post-inflationary density parameter Q is plotted for the scalar field (solid 
line) radiation (dashed line) and cold dark matter (dotted line) in the quintessential- 
inflationary model decribed by 1341 with p = 0.2. Late time oscillations of the scalar 
field ensure that the mean equation of state turns negative (w^) ~ —2/3, giving rise 
to the current epoch of cosmic acceleration with a{t) oc and present day values 
fi(i^ ~ 0.7, fiom — 0.3. From Sahni, Sami and Souradeep |rBlJj . 



The rationale for the f, R term is that quantum effects associated with 

J brane ^ 

matter fields are likely to give rise to such a term in the Einstein action as 
discussed by Sakharov in his development of induced gravity .176i. 
The resulting Hubble parameter in the DDG braneworld is 

where = rri^ /M^ is a new length scale determined by the four dimensional 
Planck mass m and and the five dimensional Planck mass M respectively. An 
important property of this model is that the acceleration of the universe is not 
caused by the presence of any 'dark energy'. Instead, since gravity becomes 
five dimensional on length scales R > Ic = 2Hf^^{l — fini)~^, one finds that 
the expansion of the universe is modified during late times instead of early 
times as in the RS model. 

A more general class of braneworld models which includes RS cosmology 
and the DDG brane as subclasses was developed in |Snil89| and is described 
by the action 



S = M^ f in- 2ylb) + f {m^R -2a)+ f i^^atter • (43) 

J bulk J branc J brane 
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For (T ~ Ay^ — (|43|l reduces to the action describing the DDG model, whereas 
for TO = it describes the Randall-Sundrum model. 

As demonstrated by Sahni and Shtanov |17(Jj the braneworld which follows 
from the action (|43|l describes an accelerating universe at late times with the 
Hubble parameter 



Hi 



2^i^n,,,{i+zY + + + , (44) 



where 



' " llHl ' " 3TO2ij2 ' - 3™2^2 . ^^A, - • I45j 

(The =F signs refer to the two different ways in which the brane can be embed- 
ded in the bulk, both signs give rise to interesting cosmology |170| .') As in the 
DDG model Ic ^ H^^ if M ~ 100 MeV. On short length scales r < and 
at early times, one recovers general relativity, whereas on large length scales 
r ^ Ic and at late times brane-related effects begin to play an important role. 
Indeed by setting M = (12; = 0) (gH) reduces to the LCDM model 



I2„,(l + z)3 + |?, , (46) 



'0 



whereas for a — — (|44|) reduces to the DDG braneworld. An important 
feature of the braneworld 144|l is that it can lead to an effective equation of 
state of dark energy WcS < — 1. This is easy to see from the expression for the 
current value of the effective equation of state jl7U| 



3(1 -12m) i-n„,^nn, + Oa + nf + nAj ^ ' 

we find that wq < —\ when we take the lower sign in H47|l , which corresponds 
to choosing one of two possible embeddings of this braneworld in the higher 
dimensional bulk. (The second choice of embedding gives > —1.) 

It is also possible, in this model, for the acceleration of the universe to 
be a transient phenomenon which ends once the universe returns to matter 
dominated expansion after the current accelerating phase. As discussed in 
| 17U| such a braneworld will not have an event horizon and may therefore 
help in reconciling an accelerating universe with the demands of string/M- 
theory. Other possibilities of obtaining dark energy from extra dimensions 
have been discussed in [7| [HI EUnOI Cni 1113 ElIISHI • The possibility that 
DE could arise due to modifications of gravitational physics has also been 
examined in [TUSl HOI HSIISSI IT^ IHSl [T^ . 
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2.5 Chaplygin gas 

A completely different route to dark energy is provided by the Chaplygin gas 
which obeys the equation of state 

Pc = -A/pc . (48) 

The conservation equation dE = —pdV ^ d{pa^) — —pd{a^) immediately 
gives 



p, = = ^A + B{\ + zf , (49) 

where i? is a constant of integration. Thus the Chaplygin gas behaves like 
pressureless dust at early times and like a cosmological constant during very 
late times. 

The Hubble parameter for a universe containing cold dark matter and the 
Chaplygin gas is given by 



H{z) = ffo 



I2™(i + z)3 + :^^| + (i + z)6 



1/2 



(50) 



where k = Pam/^/B and it is easy to see from H5()|l that 

K = ^^^^(z oo) . (51) 

Pc 

Thus, K defines the ratio between CDM and the Chaplygin gas energy densities 
at the commencement of the matter-dominated stage. It is easy to show that 

It is interesting that the Chaplygin gas can be derived from an underlying 
Lagrangian in two distinct ways: 

• One can derive it from a quintessence Lagrangian H30|l with the potential 

m 



y(0) = ^(cosh30+^) . (53) 

• The Chaplygin gas can also be derived from the Born-Infeld form of the 
Lagrangian density 



£ = -VoVT^J^ , (54) 
where cj)^^ = dcfi/dx^ . For time-like 0.^ one can define a four velocity 

= / , (55) 
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this leads to the standard form for the hydrodynamical energy-momentum 
tensor 

T^,l, ^ {p + p)Uf,u„ - pg^,^ , (56) 

where [7S] 

^" P = -VbVl-0,M'/"^ - (57) 



i.e. we have recovered ()48(l with A ^ Vq. 

The fact that the properties of the Chaplygin gas interpolate between 
those of CDM and a yl-term led to the hope that the CG might provide a 
conceptual framework for a unified model of dark matter and dark energy. It 
should however be noted that in contrast to CDM and baryons, the sound 
velocity in the Chaplygin gas Vc = y/dpc/dpc = VA/pc quickly grows cx 
during the matter-dominated regime and becomes of the order of the velocity 
of light at present (it approaches light velocity asymptotically in the distant 
future ). Thus, when one examines classical inhomogeneities, the properties 
of the Chaplygin gas during the matter-dominated epoch appear to be rather 
unusual and resemble those of hot dark matter rather than CDM, despite the 
fact that the Chaplygin gas formally carries negative pressure [2]- 

A 'generalized Chaplygin gas' has also been proposed for which p cx —1/p". 
The equation of state in this case is 

which interpolates between w = at early times (a ^ 1) and w = —1 at late 
times (a ^ 1); Wo is the current equation of state at a = 1. (The constant 
a regulates the transition time in the equation of state.) WMAP, supernovae 
and large scale sructure data have all been used to test Chaplygin gas models; 

see EziEniizHiEiiniEainsiinHiEniEsi 



2.6 Is Dark Energy a Phantom ? 

In an influential paper Caldwell noticed that a very good fit to the 
supernova-derived luminosity distance was provided by dark energy which 
violated the weak energy condition so that w < —1. He called this Phantom 
dark energy.^ Indeed, a study of high-z Sn finds that the DE equation of 
state has a 99% probability of being < —1 if no priors are placed on f2„i ' 
When these Sn results are combined with CMB and 2dFGRS the 95% confi- 
dence limits on an unevolving equation of state are —1.61 < w < —0.78 |99| . 
which is consistent with estimates made by other groups |193[ I201| . 

A universe filled with Phantom energy has some interesting but bizarre 
properties. 

* Phantom takes its name from Part I of the Star Wars movie series - the Phantom 
Menace. 
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If icq marks the epoch when the densities in matter and phantom energy 
are equal then the expansion factor of a universe dominated by phantom 
energy grows as 



a{t) ~ a (tcq) 



(l + w) 



t 



w 



cq 



2/3(l+u;) 



, W < -1 , 



and therefore diverges in a finite amount of cosmic time 



ait) 



oo as t 



w 



1 



(59) 



(60) 



By substituting w < —1 into (|24|l we immediately find that the Hubble 
parameter also diverges as t — s- teR, implying that an infinitely rapid 
expansion rate for the universe has been reached in a finite time. The 
divergence of the Hubble parameter is associated with the divergence of 
phantom density which grows without bound 



Pit) 



[1 + w) w 

tcq 



(61) 



and reaches a singular value in a finite interval of time p{t) — > oo, t — > isR- 
Thus a universe dominated by Phantom energy culminates in a future 
curvature singularity ('Big Rip') at which the notion of a classical space- 
time breaks down. (See also [Ml El Ei [HTl El 112 ED El IMI EHI HOI 
Ej.) 

The ultra-negative phantom equation of state suggests that the effective 
velocity of sound in the medium v — y^\dp/dp\ can become larger than 
the velocity of light in this model. 

Although a dynamical model of phantom energy can be constructed with 
the 'wrong' sign of the kinetic term, see (p?T|l . such models are plagued with 
instabilities at the quantum level ^Hl which makes their existence suspect. 
It should be pointed out that phantom is not the only way to get w < —1. 
A model with similar properties but sharing none of phantom's pathologies 
is the braneworld model of |170[ ^ , which has WcS < — 1 today but does 
not run into a 'Big Rip' in the future. 



2.7 Reconstructing Dark Energy and the Statefinder diagnostic 

In view of the considerable number of dark energy models suggested in the 
literature, it becomes meaningful to ask whether we can reconstruct the prop- 
erties of DE from observations in a model independent manner. This indeed 
may be possible if one notices that the the Hubble parameter is related to the 
luminosity distance |195[I17^ 



(62) 
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and that, in the case of quintessence, the scalar field potential as well as its 
equation of state can be directly expressed in terms of the Hubble parameter 
and its derivative [1951 1175| 

(63) 

(64) 
(65) 
(66) 

Both the quintessence potential V{(l>) as well as the equation of state w^{z) 
may therefore be reconstructed provided the luminosity distance dL{z) is 
known to reasonable accuracy from observations. 



SttG /d^y _ 2 dlnH n^x 
p _ {2x/3)dliiH/dx-l 
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Fig. 9. The relative difference between the Hubble parameter reconstructed from 
Sn data and the LCDM value is shown as a function of redshift. Sn data from Tonry 
et al (2003) were used for the reconstruction. The best-fit is represented by the thick 
solid line assuming fim ~ 0.3. The light (dark) grey contours represents the la {2a) 
confidence levels around the best-fit. The dashed horizontal line shows LCDM. From 
Alam, Sahni, Saini and Starobinsky 



In practice it is useful to have an ansatz for either one of three cosmological 
quantities: dL{z), H{z) or w(z), which can then be used for cosmological 
reconstruction |175[ll3lllll2[|^TT] . Popular fitting functions discussed in the 
literature include: 

(i) An ansatz for the dark energy |173j 
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N 

Pde{x) ^^^Aix\ x^l + z. (67) 

i=0 

(ii) Fitting functions to the dark energy equation of state |212[I1U7| : 

N 

w{z) = ^ WiZ"- , 
i=0 

w{z) =wo + — — . (68) 

_L ~T~ Z 

The fitting parameters Wi, Ai are obtained by matching to observations. In 
practice the first few terms in either series (|67ll , H68|) is sufficient since the cur- 
rent Sn data is quite noisy; see |46l 12121 [?)4 77, 113 , 2 for a discussion of these 
issues. An example of cosmological reconstruction of the Hubble parameter 
from Sn data is shown in figure iPJ; see also j2()8l 1134] . 

The Sn inventory is increasing dramatically every year and so are increas- 
ingly precise measurements of galaxy clustering and the CMB. To keep pace 
with the better quality observational data which will soon become available 
and the increasing sophistication of theoretical modelling, a new diagnostic of 
DE called 'Statefinder' was introduced in |173j . 

The statefinder probes the expansion dynamics of the universe through 
higher derivatives of the expansion factor a and is a natural companion to 
the deceleration parameter which depends upon d (|20|l . The statefinder pair 
{r, s} is defined as follows: 

a ^ 9w ^ , , 3 ^ w , , 

.^^ = l + -/2x(l + H-2^x^, (69) 

r — 1 1 w , , 

^"3(^3T7^-^ + "-3^- 

Inclusion of the statefinder pair {r, s}, increases the number of cosmological 
parameters to four^: H, q r, s. The Statefinder is a 'geometrical' diagnostic 
in the sense that it depends upon the expansion factor and hence upon the 
metric describing space-time. An important property of the Statefinder is that 
spatially flat LCDM corresponds to the fixed point 



{r,s} 



= {1,0}. (71) 

LCDM 



Departure of a given DE model from this fixed point provides a good way 
of establishing the 'distance' of this model from LCDM 0. As demonstrated 
in |173[ l^ r751l216j the Statefinder can successfully differentiate between a wide 
variety of DE models including the cosmological constant, quintessence, the 
Chaplygin gas, braneworld models and interacting DE models; an example is 
provided in figure ITUl 

® r has also been called 'cosmic jerk' in |2U7| 
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Fig. 10. The time evolution of tlie statefinder pair {r, s} for quintessence models 
and the Chaplygin gas. Solid lines to the right of LCDM represent tracker potentials 
V = Vo/(j>°', while those to the left correspond to the Chaplygin gas. Dot-dashed 
lines represent DE with a constant equation of state w. Tracker models tend to 
approach the LCDM fixed point {r — l,s — 0) from the right at t oo, whereas 
the Chaplygin gas approaches LCDM from the left. For Chaplygin gas k is the ratio 
between matter density and the density of the Chaplygin gas at early times. The 
dashed curve in the lower right is the envelope of all quintessence models, while the 
dashed curve in the upper left is the envelope of Chaplygin gas models (the latter 
is described by «: = ^2m/l — f^m )- The region outside the dashed curves is forbidden 
for both classes of dark energy models. The ability of the Statefinder to differentiate 
between dark energy models is clearly demonstrated. From Alam, Sahni, Saini and 
Starobinsky 

2.8 Big Rip, Big Crunch or Big Horizon ? The fate of the 
Universe in dark energy models 

The nature of dark energy affects the future of our Universe in a very signif- 
icant way. If DE is simply the cosmological constant, then the universe will 
accelerate for ever. Of great importance is the fact that an accelerating LCDM 
universe developes an event horizon similar to the one surrounding a black 
hole |196| . Consider an event (ri, ti) which we wish to observe at our location 
at r = 0. Setting ds"^ = we get 




(72) 
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Any event in the universe will one day be observed by us if the integral in the 
RHS of l|72|l diverges as t — s- oo. For power law expansion this clearly implies 
a cx i^, p < 1, i.e. a decelerating universe. In an accelerating universe exactly 
the opposite is true, the integral in the RHS converges signalling the presence 
of an event horizon. In this case our civilization will receive signals only from 
those events which satisfy |172j 

dr f°° cdt' . 

< / 777n- (73) 



VI - - Jt, a{t') " 

For de Sitter-like expansion a — aiexp H{t — ti)^ H = \J /1/3, we get r\ = 
c/aiH, so that the proper distance to the event horizon is Rh — airi = c/ H . 
In LCDM cosmology, 



H = H{t^oo)^ ./Aji - i/oVl - , (74) 
and the proper distance to the horizon is 

Rh = ^ J—r^ - 3.67/1-1 ^ ^^^^ 

if — 1/3. Thus our observable universe will progressively shrink as as- 
trophysical bodies which are not gravitationally bound to the local group 
get pushed to distances beyond Rh- (More generally, horizons exist in a 
universe which begins to perpetually accelerate after a given point of time 
| 81ll68[[l71| . To this category belong models of dark energy with equation of 
state —l<w< —1/3, as well as 'runaway scalar fields' |198| which satisfy 
V, V, V" and V'/V, V" /V ^ as ^ oo.) 

The presence of an event horizon implies that, at any given moment of 
time there is a 'sphere of influence' around our civilization. This sphere 
has an associated redshift z^, and a celestial body having z > will be 
unreachable by any signal emitted by our civilization now or in the future; 
zh — 1.8 in LCDM cosmology with Sl^ ~ 2J7m — 2/3. Thus all celestial 
bodies with 2; > 1.8 lie beyond our event horizon and there is no possibility 
of causal contact with any of them. 

Interestingly, an N-body simulation tracking the future of an LCDM uni- 
verse has shown that ~ 100 billion years from now the observable universe will 
consist of only a single massive galaxy within our event horizon - the merger 
product of the Milky Way and Andromeda galaxies |130| . Furthermore, since 
the growth of large scale structure freezes in an accelerating universe, the mass 
distribution of bound objects will cease to evolve after about 30 billion years. 

This somewhat gloomy future scenario is not absolutely essential and can 
be avoided if the currently observed acceleration of the universe is a transient 
phenomenon.^ Just such a possibility exists in a class of braneworld models 

^ Accelerating cosmologies without event horizons are important in a different con- 
text. Since the conventional S-matrix approach may not work in a universe with 
an event horizon, models with horizons may pose a serious challenge to a funda- 
mental theory of interactions such as string/M-theory. 
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| 170) in which the current accelerating phase is succeeded by a decelerating 
matter dominated regime. Quintessence potentials can also have this property, 
as discussed in [2] . An interesting class of transiently accelerating DE models 
is constructed around a scalar field potential which decays with time and 
becomes negative at late times jTH |S 11391 EHI El El- An example is V = 
Vq cos 4>/f which describes axionic quintessence jSlEHlEniEl- Such a universe 
recollapses in the future when H{to + AT) = 0, and contracts thereafter 
towards a 'Big Crunch' singularity. Supernova observations indicate that, for 
typical decaying potentials, the universe will not collapse for at least AT ~ 20 
Gyrs IS]. 

DE models have also been proposed which encounter a 'quiescent singu- 
larity' while expanding. At the 'quiescent singularity' the second derivative of 
the expansion factor diverges while its first derivative remains finite ^190. .TH] 
(i.e. a^— oo, a~ constant). In such models the expansion of the universe 
'brakes' to a virtual standstill as the universe approaches the singular regime 
at which invariants of the space-time metric diverge {RikimR^'^''™ oo) while, 
curiously, the Hubble parameter and the energy density remain finite. Cosmo- 
logical consequences of models which encounter a future quiescent singularity 
(or a 'Big Break' [JHI) have been briefly discussed in [Tl 11901 17!^ but need to 
be examined in more detail. 

Finally, as discussed in section (TBI Phantom models with w < —1 ex- 
pand towards a Big Rip, at which the density and all curvature invariants 
become infinite. As in the case of the Big Crunch singularity, the Big Rip 
will occur only in the very distant future (if it occurs at all). For instance, if 
w = constant > —1.5, Hq — 70km/sec/Mpc and fim — 0.3, the time to the 
Big Rip exceeds 22 Gyr [TM) . 

3 Conclusions and future directions 

From the theoretical standpoint the single most important question to be 
asked of dark energy is 

Is w = — 1 ? 

Rephrased in terms of the Statefinder diagnostic the question is: 

Is a /aH^ = 1 ? 

If future observations do answer this question in the affirmative^ then, in all 
likelyhood the cosmological constant is the vacuum energy, and one will need 
to review the cosmological constant problem again, in order to fathom why 
the formally infinite quantity {Tik) is in fact so very small. 

If on the other hand, either u; ^ — 1 or if the DE density is shown to be time 
dependent, then the cosmological constant problem may need to be decoupled 
from the DE conundrum and searches for evolving DE models which produce 

* i.e. if TO = —1 is measured to satisfyingly high accuracy 
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Fig. 11. Target statistical uncertainty of the SNAP experiment is sliown overlayed 
with current results from CMB and LSS observations. From Aldering 



Pde — lO^^'^GeV"' without exacerbating 'cosmic coincidence' will need to be 
examined deeply in the light of developments both in high energy physics and 
in gravitation theory (superstring/M-theory, extra dimensions etc.). In either 
case the key to determining the properties of DE to great precision clearly lies 
with ongoing and future astrophysical experiments and observations. 

Since the original discovery of an accelerating universe |152l 11531 1159| the 
Sn data base has grown considerably and data pertaining to 200 type la 
supernovae are avahable in the hterature [221 IMl CHI HEDl • Although system- 
atic effects such as luminosity evolution, dimming by intervening extragalactic 
material (alternatively brightening due to gravitational lensing) continue to 
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be a cause of some concern - recall that a luminosity evolution of ~ 25% over 
a lookback time of ~ 5 Gyr is sufficient to nullify the DE hypothesis |158j - 
it is reassuring that recent observations of CMB anisotropics and estimates of 
galaxy clustering in the 2dF and SDSS surveys, make a strong and indepen- 
dent case for dark energy |193[l200ll201 |. Indeed, a joint analysis of CMB data 
from WMAP + HST Key Project determination of Hq imply wde < —0.5 at 
the 95% confidence level |19Hj . 

It is of paramount importance that Sn observations continue to be supple- 
mented by other investigations which are sensitive to the geometry of space 
and can be used as independent tests of the DE hypothesis. The volume- 
redshift test, Sunyaev-Zeldovich surveys, the Alcock-Paczynski test, the an- 
gular size-redshift test and gravitational lensing have all been suggested as 
possible probes of dark energy, and will doubtless enrich the theory vs obser- 
vations debate in the near future. In addition, the proposed SNAP satellite 
which aims to measure light curves of ~ 2000 supernovae within a single year 
| 219| . should provide a big step forward in our understanding of type la su- 
pernovae and help determine the cosmological parameters to great precision, 
as shown in figure ITTl 
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